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Abstract 

The spiner fermulation of magnetogas dynamics that was developed by 
the author in a previous paper is applied t© the problem of a sound 
wave in an electrically neutral gas. This simple, one-dimensional ^problem 
serves t# illustrate the essential features ef the formalism. The solution 
is completely relativistic and, fer sound waves of macroscopic wavelengths, 
satisfies Euler’s equation, the continuity equation, and the adiabatic 
condition. For wavelengths of the order of atomic dimensions, spin-dependent 
terms, which are completely negligible for macroscopic waves, become 
important and drastically alter the form of the solution. Thus the spinor 
solution of the sound wave problem has the desirable feature that it 
automatically breaks down at the point where classical mechanics breaks 
down, namely when the characteristic length of the problem 
atomic dimensions. 

I ntroduction 

This paper is a supplement to the author’s earlier paper "Spinor 
Formulation of Magnetogas Dynamics" (Goddard Space Flight Center X-6i|.0-6i|-55) . 
In order to illustrate the formalism developed in that paper (which shall 
henceforth be referred to as "I"), the spinor formulation of fluid dynamics 
is applied to the one-dimensional problem of a sound wave in an electrically 
neutral perfect gas. 

The solution is carried out by means of a perturbation on a simple 
zero-order solution which corresponds to a fluid of constant density 
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moving in the positive z direction with constant velocity. By choosing 
the particle spins, whose role in any macroscopic problem is always 
insignificant, to be aligned in the z direction, it is possible to reduce 
two of the four spinor components to zero, thereby reducing the calcula- 
tional burden. The fluid enthalpy, which is regarded as the driving func- 
tion of the problem, is next assigned a sinusoidal time-independent vari- 
ation in the z direction, and the spinor equations are solved to find the 
spinor functions consistent with this form of variation in the enthalpy. 

The solution is carried out using the first-order perturbation approxima- 
tion in which the perturbation parameter is the ratio of the maximum change 
in the specific enthalpy to the particle rest-mass. For practical problems, 
this ratio is always very small. 

Once the spinor equations have been solved for a sinusoidal variation 
in enthalpy, spinor relations derived in I are applied to calculate the 
fluid flux density, which turns out to be constant, and the particle density, 
which turns out to have a sinusoidal variation of the same wavelength as 
the assigned variation in the enthalpy. It is shown that the solution 
maintains conservation of particle energy, which is just the condition 
required by Euler's equation. That this is so can be seen from the fact 
that for this problem eq. ( 2 - 37 ) of I (Euler's equation), reduces to 

0=V( V^mc 2 - + rn'fc.) ( 1 ) 

which states that the sum of the kinetic and thermal energies per particle 


must remain constant 



The spinor equations, like Euler's equation, must be supplemented 
by the adiabatic condition (eq. (2-31) of I): 

ff-R. = Or-i) & 

where and S & are the changes in the density and enthalpy respectively 
as one moves along the wave, which appears to be stationary because the 
fluid is streaming in the+z direction with exactly the same speed as that 
with which the wave is propagating in the -z direction. It is shown 
that (2) yields a condition on the fluid velocity which is just the usual 
expression for the speed of sound in terms of the absolute temperature of 
the gas. 

Finally, it is shown that, although the effects of particle spin are 
completely negligible for macroscopic wavelengths, they become important 
when the wavelengths become of atomic dimensions, with the result that the 
classical solution is no longer valid. 

Necessary Spinor Relations 

The spinor relations that will be needed are recapitulated below. 

The numbers to the left are the formula references in I . Since in the 
present problem two of the spinor components ( and Z ) are zero, the 
spinor equations take the following simplified form: 


■la,b) 

i(± 



li 

1 

i 


if— 



- )r 


L\ c. 

TF ~ 

Ti" ) 
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where the overhead bar indicates complex conjugation and 


h 


where 

(8-ic) tri - m(M- &S ? ) (5) 

c 4 

where is the variable part of the enthalpy, the constant part having 
been absorbed into the particle rest-mass m. 

Once the spinors are known, the flux density U -vector is calcu- 

lated from the following relations: 


(B-15) 


pu°=^[is'r+i£T+i¥'r+ 




fy' = re tRlr'l 2- + %' z 1 ' ] = o 

fu 1 --yr5:4'|YF + x'x T ] - o 


( 6 ) 


•J r tif , r-ifT+ix , r-ixTj =j=[i?V'-i*v] 

x I 

where the simplifications that result when ? - X — O have been carried 
out. The invariant density is found from the relation 


(8-9) C = t(ey J )<;etA i )]' /, -= = \fz U'x 1 ! (7) 


The particle 3-velocity /v * is given by the following well-known relation: 

= C (^«) J = 1,2-, 3 

Finally, the particle li -momentum -jp) is found as follows: 

IJ . )] 


(8) 


• (C-Sb) 

( c— 79 ) 


( 9 ) 
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Zero-Order Solution 


The zero-order function is taken t© have 


the following form: 


n = o 

X*. * O 
Xo = B© e 


l(Et-P*)/fc 


( 10 ) 


where A q , B q , E, and P are constants. (10) is a solution of (3) for £■&= O 
if the following relations are satisfied: 


B© = 

A© — 


E-Pc 

me 1 


Ac 


E-hPc 
no c** 


Bo 


(lla) 

(lib) 


These can be satisfied only if 

E* = (me 1 ) 1 +- (Pc)* - (12) 

which is the familiar relativistic relation between the energy E and the 
momentum P of a particle having rest-mass m. 

A 

It will be convenient to work with the dimensionless quantities E and 

A 

P defined as follows: 




A 

E 

£ 

me** 

_ 1 
VT~ <y/c.) x 

(13a) 



A 

P 

_ _P_ = 
rvic 

(V/c) 

]f\ - (v/c.)* 

(13b) 

Thus (12) assumes 

the 

following 

form: 





V- 

A i 

= 1 +- P 

L 

(lU) 





B c 


For convenience we 

introduce the 

ratio R - 

! A* : 



R - 


A 

-P = ~ 

1 

A 

(15) 



6 


The following identities will be useful: 


|+R l = IKE 


( 16 a) 

\- R l = iRP 


(16b) 

Substituting (10) into (7) and (6), we find 

e. = fi R 


(17) 

^(i +R ’') = 

A 

foE 

(18a) 

(eu 3 ) c = = 


(18b) 

Substitution of (18) into (8) yields 

3 P P 

^ " C E " (P/c*) 

= V 

(19) 


This equation is taken as the definition of the constant velocity V, which 
will later be identified with the velocity of sound. 


We note that the nonrelativistic case is defined by the following conditions: 


A 


P = 


_P_ 

me 


« s 


A 


E = 



N.R. case 


( 20 ) 


In this case 


n « V ~ £ 


N.R. case 


( 21 ) 


Solution in the Presence of a Sound Wave 

In order to generate a sound wave, we regard the change in specific 

enthalpy as the driving function, and assign it the following functional form: 

& =■ A cos k* (22) 

where k is the wave number of the sound wave. (For an ordinary sound wave 

—2 —l 

in hydrogen at S.T.P. k ^ 10 cm .) £ "fo appears in the spinor equations 

as a contribution to m: 


m = m (i ~ ) = m[\ -t- (A)cos ki] 


(23> 


The dimensionless constant ^ plays the role of the perturbation para- 
meter. For a strong sound wave at S.T.P., for which the variation in 
pressure produced by the wave is of the order lO^dyne/cm^, the constant ■— 

-15 

is ef order 10 . 

Because we wish te satisfy the spinor equations (3) only to first order 
in A , we can approximate a term like by $ 0 (S£). 

We seek a solution of (3) having the form 

= A^l + <x C05 HE 3- k z] e*^ Et ? 



t r -i -t(Et-PE)A. 

X - B 0 [l 4-/Q cos kH: - t(§ + 7?) sin ks] e * 

where o< , ^3 , § , and ^ are constants. 

The constant ^ has "the effect of producing phase changes in § and % 
that are equal in magnitude and opposite in sign. When (2U) is substituted 
into (6) and (7), % dees not appear to first order in A , It is only 

in (9) that it makes itself felt, and even there the effect is negligible 
for macroscopic waves. In fact, as an alternative procedure, we could omit 
£ entirely and compensate for the omission by replacing Pz in the phase 
by an arbitrary function of z. We shall later see what this function would 
have to be. 

The constant ^ is a measure of the degree to which the constraint 

is violated. It was pointed out in Section 8 
of I that, although in principle this constraint should be imposed, the 
effect of neglecting it in any macroscopic problem is completely negligible. 
We shall see that this is true in the case of a sound wave. 
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It is to be noted that the two constants ^ and , being coefficients 

I X 

of the sin kz term, specify the part of f and X that is out of phase with 
sn, whose variation is given by A cos kz. We shall see that these constants 
have no influence on macrescopic solutions. Only the in-phase constants c< 
and ^ are important for such solutions. 

The spinor specification given in ( 2k ) is a solution of (3) (to first 

4. ^ ^ 4 f rt«rie+Qn+.e ViMr. + V\o f.11 mH tut lrnl 11 .C • 


where 


expressions! 


o _ » /A \ T 1 tP(E-»-?) 1 
P ~ X 

(25a) 

(25b) 

S - _L ( A_\ pp ‘KVO* 1 

> * u ; )PL p--(^J 

(25c) 

^ ^ E (A/c*) 

(25d) 

V e M » A f 

»= k X 

(25e) 

(2^), and the relations (6)-(8), we arrive at the 

following 

f = Vi A^R^I + Ox-M) »> k*] 


«l 

fo ^ 1 + C c«.) ^ pt. _ (9/-^* ] C ° 3 ^ 

(26) 


- o t T| 4- (A/C ^ C0S 

^^(»-R x )[\-t-Z(^-R a >3)cos kzjj 

-e.p 


(27) 


( 28 ) 


( 29 ) 
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u°= 


e = 4- (\- R l ) (*-/3) cos k*] 

= g eos **} 

From (23) and (30) we find 

t (to r~ y /x)X j cos k Z = E 

c *• ? ~(®A) J 

All of the above results are, ef course, valid only to first order in ~ 
For all macroscopic sound waves ©« | . In fact, for an ordinary 
sound wave in hydrogen k'V'lO 2 cm _1 and *v lO'^cm, so in this case 
& SV lO -1 ^. For comparison, in the 


* P 

P* 


m a 


same case 


£ - IO-' 


Thus (?/%) is completely negligible compared with P 2 . Using this fact 
in (26) -(31), we have 


e = e„ jj 




(30) 


(31) 


(32) 


(33) 


+ cp/^] 

(w 5 = e.c^) 


(3U) 

(35) 

(36) 

(37) 


1* O 

me H — El 


(38) 


These expressions are, of course, completely relativistic 


10 


From (27) and (28) we see that 


d j* ) - c •+" = O 


'J" ” ' C. Int ~ y ^9) 

Thus the continuity condition is always fulfilled, even when ©'is not negligible. 

(38) is obviously a solution to (l), which is just Euler's equation for the 

case of a sound wave. In (31), however, we find another terra appearing in 

the energy equation. For a strong sound wave in hydrogen at S.T.P. for 

-? -1 - t K*- 

which k^lO cm , this term has a magnitude of the order 10 erg as 

-3 o 

compared with a magnitude of order 10 erg for the term me u . This 
small extra energy contribution results from the presence of particle spin. 

For macroscopic sound waves it is completely negligible, but for extremely 
short wavelengths it becomes important. 

We shall return to th question of short wavelengths, but first we 
shall impose the adiabatic condition (2) in order to arrive at the familiar 
expression for the speed of sound. We note that^to first order in , 

it is permissible to write (2) as follows: 


e. si = (v-'K.Je = (*-0-*,% se (to) 

From (33) we have 

(1.1) 

V 

Substituting (Ul) into (I 4 . 0 ) we have 

V = fn-iKpTo 0 , 2 ) 

which is the well-known expression for the speed of sound in a perfect gas^ 
where is the constant-pressure specific heat per unit mass and is 
the ratio of specific heats. 

To complete the solution, we apply (?) to arrive at the following 
expression for particle momentum: 
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‘t 1 '* - P - ti k % cos k£ 


(l;3a) 



P v -(<5/0* 


0i3b) 


For the macroscopic case, the second term on the right side of (U3b) is 
completely negligible, and so and stand in just the linear relation 
to each other that we would expect. 

It is obvious that (l*3a) could be rewritten as follows: 

= ji ~ ^ 5 5 ' v ' k *l (Ui) 


This indicates that if we had replaced Pz in the phase factor in (2li) by 
Pz-fc|sin kz where $ has the value given in (25c), then it would have 
been unnecessary to include % in the square brackets multiplying A and 
B # in (2U). 

Spin-Velocity Resonance 

_ ^ 

The expressions (25) have the interesting feature that for < %* P or 

P = K U5) 

the denominators vanish, so that even for A - O it would be possible to 
have a wave of finite magnitude. The physical explanation for this is 
that for very short wavelengths, the spin-dependent forces that are negHk 
gible for macroscopic waves become large enough to play the role that is 
played by the pressure in a macroscopic sound wave. This phenomenon could 
be called spin-velocity resonance. It occurs only for wavelengths of the 
order of atomic dimensions, for which classical, mechanics loses its vali dity. 

For example, for waves in hydrogen at S.T.P., resonance would occur at a 
wavelength of the order 10 "^ cm. 

When resonance occurs, the solution given above is no longer valid 
(even if the ratio *-r — — - — were held constant) because the perturba- 
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tion procedure is no longer valid. Moreover, before any physical reality 
could be claimed for such a resonance, it would be necessary to find a 
solution that satisfied the constraint J- (?* Z") ~O f because otherwise 
the forces arising from the neglect of this constraint, which are completely 
spurious and without physical meaning, would play an important role. 

The occurrence of resonance does, however, represent a breakdown of 
the classical solution and so serves to answer the question "Why, if the 
constant 1i plays n© role in macroscopic problems, must it be chosen to 
have the value 1.0!?X10“ 'erg-sec if our only interest in the spinor alterna- 
tive to Euler* s equation is to solve macroscopic problems?” The answer 
is that, in order to give the correct breakdown point, namely the point 
at which the characteristic length of the problem becomes of atomic 
dimensions, the value assigned to in the classical spinor theory must 
be at least approximately equal t© the value given above. 
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Note Concerning Speed of Sound 

It should be noted that^in the expression (i|2) for the speed of sound 
is a function of the temperature T 0 . This is most directly seen from the 
expression for in terms of Boltzmann's constant k , particle mass m, 
and number of degrees of freedom f: 


But the mass m contains a contribution from the enthalpy , and so is 
a function of T 0 . Now let /Cp be the specific heat referred to the particle 
mass m Q that does not contain the enthalpy contribution: 

= -s o * i) 

mt l = m„c x +- vr > e ^ 0 - m e> (c 7 ' T 0 ) 


0 * 6 ) 


Since 


we have 


r - rn ° ^ - gX/0 P 

A-p - /°p ^ rn ) ~~ r » +/ o T 


C*+/C*T e 


which when substituted into (1*5) yields 
According to (1*9) 


A 


0*7) 

m 

0 * 9 ) 

( 5 0) 

( 51 ) 


lin, V= cVOM) = c\[f 

° ->r 

This limit cannot, however, be taken very seriously since, as Synge points out 
on p.58 of his book. on the relativistic gas, equation (2) for the adiabatic condition 
is not valid at temperatures for which the thermal energy is comparable 
with the rest energy. The limiting speed of sound for a monatomic gas 
(f = 3) that Synge gives (p.77) is 

V, = ( 52 ) 

*»•*» vr 


J.L. Synge, The Relativistic Gas (Inters cience Publishers, New York, 1957) 


